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The electronic nematic order characterized by broken rotational symmetry has been suggested 
to play an important role in the phase diagram of the high temperature cuprates. We study the 
interplay between the electronic nematic order and a spin density wave order in the presence of a 
magnetic field. We show that a cooperation of the nematicity and the magnetic field induces a finite 
coupling between the spin density wave and spin-triplet staggered flux orders. As a consequence 
of such a coupling, the magnon gap decreases as the magnetic field increases, and it eventually 
condenses beyond a critical magnetic field leading to a field-induced spin density wave order. Both 
commensurate and incommensurate orders are studied, and the experimental implications of our 
findings are discussed. 
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Introduction — Typically, conduction electrons in 
strongly correlated materials are either localized due to 
interactions between them, or conduct with a uniform 
and isotropic distribution. Can electrons in metals as- 
semble themselves and exhibit novel patterns seen in liq- 
uid crystals? If so, what is the mechanism behind such a 
self-organizing pattern, and what are the effects of such 
a metal on nearby phases? The electronic nematic phase 
characterized by an anisotropic conduction has been pro- 
posed as one such state. It was suggested that it plays 
a relevant role in determining the superconducting tran- 
sition temperature in the cuprates which motivated 
several studies on the interplay between the nematic and 
d-wave superconducting (dSC) states. @,0| An evidence 
of its existence in the cuprates was found by neutron scat- 
tering on YBa 2 Cu 3 6 .4 5 (YBCO).Q On the other hand, 
the interplay between the nematic and the antiferromag- 
net, another nearby phase of the cuprates, has not been 
addressed. 

In this paper, we study the relationship between the 
nematic and spin density wave (SDW) orders. It is triv- 
ial to see that a direct coupling term between the two is 
not allowed in the free energy, since they break distinctly 
different symmetries - the nematic order breaks the ro- 
tational symmetry between x and y directions, while the 
AF order breaks the translational, spin-rotational, and 
time-reversal symmetries. Thus, a naive conclusion is 
that there is no strong influence between the two. We 
show that the situation can be dramatically different 
when the time reversal symmetry is broken by an ex- 
ternal perturbation such as a magnetic field. We find 
that there is a direct coupling between the SDW and 
spin-triplet staggered flux (tSF) orders inside the nematic 
phase when the magnetic field is applied. The tSF phase, 
like the (spin-singlet) staggered flux (or sometimes called 
d-density wave) breaks the translational and rotational 
symmetries due to circulating currents with an alternat- 
ing pattern. However, unlike the d-density wave state, 



it does not break the time reversal symmetry, because 
up- and down-spin circulations have opposite directions, 
which leads to the name of spin-triplet staggered flux. A 
consequence of such a coupling is that the SDW order can 
be induced by the magnetic field via the nematicity. This 
result applies to both commensurate and incommensu- 
rate orders - a finite coupling between incommensurate 
spin density wave (ISDW) and incommensurate triplet 
staggered flux (ItSF) orders occurs in the presence of a 
magnetic field leading to ISDW order beyond a critical 
magnetic field. We will discuss the experimental implica- 
tions of our result in the context of the high temperature 
cuprates. 

Nematicity and commensurate SDW and tSF orders 
- The nematic order, which breaks the 90 rotational 
symmetry between x- and y-directions, is defined as[^,0, 

a 

A/o = ^d(k)4 CT c kff , (1) 

k<x 

where d(k) = cos k x — cos k y (setting the lattice constant 
a = 1). Inside the nematic state, marked by (a) in the 
phase diagram of Fig. [TJ the quasiparticle Green's func- 
tion is written as 

G _1 (k, iuj n ) = -iu n + e k - /i, (2) 

where 

£k = — 2t(cos k x + cos k y ) + 2td(k)J\fo — 4t' cos k x cos k y , 

and [L is the chemical potential, t and t' represent the 
nearest neighbor and second nearest neighbor hoppings, 
respectively. 

As shown in Fig[T] the commensurate antiferromag- 
netic (AF) phase exists below a certain doping of x Cl and 
its order parameter with a moment along x is given by 

Sq = J2 c *o 1 <0 c *+QP> ( 3 ) 
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FIG. 1: Schematic phase diagram of SDW and nematic phases 
as a function of doping x, temperature, T and magnetic field 
B. The lines in the boxes are the dispersions of collective 
modes, (a) In the nematic state, both the SDW and tSF 
modes are gapped, (b) As the magnetic field increases, the 
magnon gap decreases, while the tSF gap gets pushed up. 
Further magnetic field makes the magnon condense leading 
to a field-induced SDW (fSDW) in the yellow region. See the 
main text for further discussions on detection the collective 
modes for incommensurate and commensurate SDW orders 
via neutron scattering and on the anisotropic intensity be- 
tween k x and k y due to the nematicity. 



where a is the Pauli matrix, i = x,y, z, and Q = (n, 7r). 
It is straightforward to check that there is no direct 
coupling between the nematic and AF orders based on 
the distinctly different quantum numbers associated with 
their order parameters. 

However, when we consider the following tSF order 
parameter, the situation changes: 
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The tSF order is characterized by circulating currents 
with an alternating pattern like the d-density wave state, 
but up- and down-spins circulate in opposite directions. 
Thus, there is no net charge current, but a finite spin 
current. It breaks the translational, x — y rotational, and 
spin rotational symmetries, but preserves the time rever- 
sal symmetry. This state was identified as a component of 
the 6 dimensional superspin which can be rotated under 
the 15 generators forming the SO(6) groupfH which in- 
cludes as a subset the SO (5) group suggested for a unified 
theory of the high Tc cupratesQ. It was also discussed 
as one of non-zero angular momentum condensate states 
in Ref. fijl. In a dSC state, both the tSF and AF modes 



are gapped while the AF gap, i.e. the magnon gap, is 
smaller than the tSF gap as shown in the box (a) in Fig. 
[TJ due to the proximity of the AF phase in the phase 
diagram. 

Effects of a magnetic field in the nematic phase - 
Since both Tq and Sq carry charge 0, spin 1, and mo- 
mentum Q, one may expect that there is a direct coupling 
between the two order parameters such as 



T = ll3 (M ,B)S l Q T^ 



(5) 



where the coupling 7^ is a function of A/o and the mag- 
netic field B to respect the discrete symmetries of the 
system. First, the tSF order breaks x — y symmetry, 
while the AF does not, thus the coupling between them 
vanishes, except inside the nematic phase. Therefore, 
7(A/o = 0, B) = 0. Second, the AF order breaks the 
time reversal symmetry, while the tSF order does not, 
which implies 7(A/o,B = 0) = 0. Therefore, the cou- 
pling between the tSF and AF orders is finite only in the 
presence of both a magnetic field and the nematic order, 
j(M ^0,B^0)^0. To linear order in Af and B, the 
form of coupling is found to be T ex A/qB • (Sq x Tq), 
i.e. the index dependence of 7^ is given by 7^ ex eijkB k . 
It is straightforward to check that such a term is allowed 
in the free energy based on the symmetry consideration 
discussed above. 

To compute 7^ in Eq. (0, it is useful to introduce 

= ( c kt ' c <- 



written as 

H- 



kt' ^iL+Qj.)- I n ^ ne basis of -0, the Hamiltonian is 
= »kPk + B)T 3 -Mk/]^k, (6) 
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Ck + Ck + Q 
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where /i k = M ~ ~ ' 2"™ = /^k+Q, e k - 2 
— ek+Q, and B = Bz. Note that the AF and tSF fluctu- 
ations couple to fermions as 
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d ( k ) T Q c U°k+Ql + h.C 
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where g\ and 172 are the interaction strengths in the AF 
and tSF channels, respectively. We have chosen x- and 
y-component of the AF and tSF fluctuations assuming 
that there is a small anisotropy such as Dzyaloshinskii- 
Moriya interaction aligning the staggered moment in the 
x — y plane. 

Then "f xy in the nematic state and in the presence of 
a magnetic field is obtained by computing the Feymann 
diagram shown in Fig. [21 



Ixy = 9192 



E 



e k + B 



[n F (e k + B - /i k ) - n F (-e k - B - /x k )] , 



which is proportional to MqB. For example, for t 
f = -0.3, and /i = -0.867, j xy = 2.36Af B gi g 2 . 



(8) 
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FIG. 2: Diagram to compute the coupling constant y xy in 
Eq. ©. 



Field-induced antiferromagnetism, collective modes, 
and neutron scattering — The free energy deep in the 
nematic state (static No), expressed in terms of the AF 
and tSF fluctuations, is given by 



T 



2m, 



[(1 +Af )(d x S Q ) 2 + (l~Af )(d y S Q ) 2 



+ ^ [(1 + N )(d x T Q ) 2 + (1 - No){d y T Q f 



2 A.St l + -A t TZ-> r S Q T Q , 



(9) 



where A s and A t represent the gap of AF and tSF, re- 
spectively, when a magnetic field is absent. 

The above free energy becomes unstable to the forma- 
tion of antiferromagnetism when the following condition 
is satisfied: 



1^ = ^^ 



(10) 



Strictly speaking, the condensed mode is not a pure AF, 
but a mixture of AF and tSF. However, if A t ^> A s , the 
mode is dominated by AF. Assuming that A s oc (x — x c ), 
the doping dependence of the critical field for the onset 
of the field-induced AF state is given by 



B c oc 



y/{x -X c )A t 

Afo 



(11) 



which determines the phase boundary between the pure 
nematic phase and the field-induced AF phase as shown 
in the dashed line in Fig. [1] 

The collective mode dispersions are shown by different 
lines in the boxes in Fig. [1] The gapped magnon de- 
noted by the thick line in the nematic state, as shown in 
Fig- Ha), can be detected by neutron scattering, and its 
intensity should be anisotropic in momentum at a given 
frequency. The anisotropic intensity of magnetic excita- 
tions was reported in YBCO, and interpreted as evidence 
of the nematic state. [H, EH While the magnon can be de- 
tected by the neutron scattering, the tSF mode does not 
couple to the magnetic field directly. Therefore, the tSF 
mode shown in the dotted line in Fig. QJa) cannot be 
detected by the neutron scattering technique. However, 
in the presence of a field, the two modes couple as we 
have shown above. The magnon gap is pushed down, 
while the tSF gap is pushed up, and the coupling also 
generates a finite but small intensity of the tSF mode 



shown by the dashed line in Fig. [lib), so the two modes 
should be in principle detectable by neutron scattering 
in a field. A further magnetic field decreases the magnon 
gap as shown in Fig. 1(c) leading a field induced AF 
order, while the tSF gap becomes larger. 

Incommensurate SDW and tSF orders — Recently it 
was reported that the intensity of inelastic neutron scat- 
tering peaks at the incommensurate wavevectors, Q = 
(tt ± S, tt) in YBCO increases when the magnetic field is 
applied. [1 21] Here we show that the field induced SDW 
can be also applied to incommensurate orders for both 
collinear and spiral cases. We will discuss how to distin- 
guish spiral and collinear orders at the end of the section. 

Let us first consider a collinear spin density wave order. 



(12) 



where we choose the spin density wave as S x (r) oc cos(Q • 
r). Taking into account broken symmetries as before, the 
following incommensurate tSF couples to ISDW linearly 
under the magnetic field. 



k 



(13) 

where J(k, Q) = 2(smfe c sin ^ — sink y sin We used 
a symmetric form of the order parameter for convenience. 
Note that the form factor of d(k, Q) has no longer a d- 
waveness when Q deviates from (tt, tt) due to the con- 
straint of the current conservation at each site. [l3l| 

Note that due to the incommensubility 5, a deviation 
from tt, there is a finite coupling between ISDW and ItSF 
under the field, even in the absence of the nematicity. 
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Mk-s) 
(14) 



However, the coupling j xy is proportional to 5 2 which 
makes j xy to be 0.022_B when 8 = 0.12-7T. On the other 
hand, when the nematicity is finite, "f xy — 2.AMoBg\g2 
similar to AF case. 

Now let us study a non-collinear or spiral spin density 
wave order. The spiral spin density wave can be written 
as 



s Q=Ek« 



"k+Q^k-Sf 



(15) 



where S<g(r) = cosQ • rx + sinQ • ry with spins lying in 
the x — y plane. A similar spiral incommensurate stag- 
gered flux can be defined as 

T Q = E d ~( k ' Q) ( c L2 T c k+§i + Ck+Q^k-Sr) , (16) 
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Similar to the spiral spin density wave order, the spiral 
staggered flux can be viewed as a pattern of incommen- 
surate staggering current where a spin quantization axis 
shifts from site to site determined by Q. Note that the 
coupling between the spiral spin density wave and the 
spiral staggered flux in the presence of a magnetic field 
is the same as that for the collinear case. 

Our results in general support a magnetic field-induced 
spin density wave including both collinear and spiral or- 
ders. Which one of these can be finally stabilized in the 
magnetic field requires a microscopic understanding of 
cuprates which is beyond the scope of our current study. 
However, we offer an experimental way to distinguish the 
collinear and spiral orders. Note that the incommensu- 
rate collinear spin density wave couples to charge den- 
sity wave order, which implies that there should be extra 
peaks at (±2(5, 0) or (0, ±26) as the spin density wave 
order sets in at (tt ± 5, tt) or (tt, tt ± 6) under the mag- 
netic field. On the other hand, the spiral spin density 
wave does not induce a charge modulation. Therefore, a 
further neutron scattering under a higher magnetic field 
would be a way to determine the precise field-induced 
order. 

Other effects of the magnetic field — The effects of 
the magnetic field on the order parameters deserves some 
discussion. In the weak coupling approach, it was shown 
that the magnetic field enhances the nematic order for 
one spin-component, but weakens the other component, 
and it shifts the phase boundary between the nematic 
and isotropic metal. |14( However, deep inside the nematic 
state, the field effect is negligible. The magnetic field 
does not directly couple to the tSF order, because the 
tSF does not break the time reversal symmetry. Thus 
the magnetic field effect on the nematic and tSF orders 
does not affect our analysis. 

The full effect of the direct interaction between mag- 
netic field and AF order will be discussed later in the 
context of high temperature cuprates, as it depends 
on circumstances. However, the effect of the Zeeman 
coupling on the quasi-particles is rather straightforward 
to compute. When the staggered moment lying in 
the x — y plane, and the magnetic field is along z, it 
changes the electronic dispersion as E-^ = Ck+ ^ k+Q ± 

(fik — £ k+Q ± 2B) 2 + 4Sq, where there is no linear B 
dependence. However, when the staggered moment and 
the magnetic field are parallel, say along z, then the dis- 
persion is E k = ^+^ ±2B ± ^(ek - e k+Q ) 2 + 4S* , 
with the linear B dependence similar to the effect of the 
magnetic field on spin singlet condensate states such as 
a charge density wave. 

Fermi surface and Quantum oscillations — The ne- 
matic phase in the absence of magnetic fields is metal- 
lic, and thus the field-induced SDW phase is expected to 
be a metal, fill] For example, as shown in Fig. [31 in the 
field-induced AF coexisting with nematic order, there are 



elongated pockets — one electron-like and one hole-like. 
It is straightfoward to generalize the Fermi surface for 
a sprial order, where the electron pocket is qualitatively 
the same. 
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FIG. 3: The topology of a Fermi surface in a field-induced 
AF via the nematic order. Note the anisotropy between the 
electron pockets near (tt, 0) and (0, tt) due to the nematic or- 
der. It is straightforward to generalize it for incommensurate 
SDW orders, where the electron pocket remains the samefl^]. 
while the hole pocket is sensitive to the incommensubility. 

In a metal with a closed Fermi surface, magnetization 
M and conductivity oscillate in 1 /B as 



M oc A n cos 



n(A eF + <t>) 
B 



(17) 



where A €F is the area of a closed Fermi surface at the 
chemical potential. Since there are two types of pocket, 
the primary periodicities (n = 1) are determined by the 
size of both the hole and electron pockets. Using the 
same parameter set of t = 1, t' = —0.3, /i = —0.867 for 
7, and setting (Sq) — 0.07 and A/o = 0.05, we obtain 
the periodicity of 540T and 900T for the electron and 
hole-pocket, respectively, as a unit cell of square lattice 
~ 3.82 x 3.89 A 2 . Due to the nematic order, the elec- 
tron pockets around (tt, 0) and (0, tt) are elongated along 
different directions, but their area is the same yielding 
a single frequency of the quantum oscillations. As dis- 
cussed above, in a Fermi liquid, the Zeeman effect gener- 
ates two different phases and <j>i (up and down spins 
have different phases) , due to the linear B dependence in 
the electronic dispersion. However, in the field-induced 
AF with a staggered moment lying in the x — y plane, 
there is no linear B dependence, implying a single phase 
of oscillation. 

It is tempting to argue that our finding is relevant to 
the quantum oscillations observed in high temperature 
YBa 2 Cu 3 6 .5 and YBa 2 Cu 4 8 .[l2l EE El H El [22| 
Recent quantum oscillations reported in YBCO around 
10 — 12% doping were striking and contrast with the 
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Fermi arc - a truncated Fermi surface - re por ted by the 
angle resolved photoemission spectroscopy. 23, 24j There 
have been a few theoretical proposals which attempt to 
explain the observed quantum oscillations, including the 
d-density wave as a pseudogap state [25j, the staggered 
flux in a vortex liquid [2||, stripe formation 27 1, and a 
field-induced antiferromagnet[28j. 

As shown in Ref. (28l |. the competition between the 
dSC and AF orders also leads to a field-induced AF 
phase, where the coexistence of AF and dSC states be- 
tween < x < 0.1 is important to induce a staggered 
moment beyond 10% doping. In our study, the AF and 
dSC coexistence is not a necessary ingredient in achieving 
a field-induced AF order. Nonetheless, the competition 
would help to lower the critical value of the magnetic field 
at which the AF order sets in. While the mechanisms of 
the field-induced AF are not shared by the two studies, 
both results favor a field-induced AF order. 

Among the commensurate orders, it is worthwhile to 
note that both the d-density wave (spin singlet stag- 
gered flux) and AF orders lead to a qualitatively similar 
Fermi surface topology originated from the broken trans- 
lational symmetry. As proposed in Ref. [29(, the ortho- 
II potential can be used to test two different scenarios, 
the d-density wave order versus AF order. It was re- 
ported that both ortho-II-free, YBa2Cu40s and ortho-II, 
YBa2Cu3C>6.5 show a primary oscillation with a period- 
icity of 660T and 540T, which supports that the primary 
oscillation is from a pocket insensitive to the ortho-II po- 
tential. Within our picture, the electron pocket is not 
affected by the ortho-II potential. On the other hand, 
the hole pockets are modified by the ortho-II potential 
as shown in Ref. [2^, and the topologies of the modified 
hole pocket are different for the AF and d-density wave 
orders. Thus, a detailed study of the magnetic break- 
down due to the ortho-II potential would be another way 
to differentiate a field-induced AF via the nematicity and 
the d-density wave order. [30| 

In summary we show that nematicity is a route to 
achieve a field-induced spin density wave order in the 
presence of a magnetic field. An induced spin density 
wave is obtained via the coupling between the spin den- 
sity wave mode and the spin triplet staggered flux mode 
which are both gapped in the nematic phase. Such a di- 
rect coupling between the two orders occurs only when 
a magnetic field is applied, and the coupling strength in- 
creases linearly to the magnetic field. A further increase 
of the magnetic field eventually condenses the magnon 
leading to a field- induced spin density wave order, while 
the gap of the triplet staggered flux mode becomes larger. 
A neutron scattering study under a high magnetic field 
will offer a precise order among possible commensurate 
and incommensurate orders discussed above. Our study 
indicates that the nematicity plays an important role 
in understanding the phenomena observed in the high 



temperature cuprates, in particular in the context of 
anisotropic magnetic excitations and quantum oscilla- 
tions observed in YBCO materials. 
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